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摘 要 : 针对 三 维 重 建 网 格 模型 经 常 出 现 异 常 缺 失 有 孔洞 的 问题 ， 提 出 一 种 缺陷 孔洞 自动 识别 与 孔洞 区 域 细节 特征 保 
持 的 曲面 修复 方法 。 首 先 对 缺失 区 域 的 上 下 文 及 轮廓 曲线 进行 异常 检测 以 判断 是 否 为 缺陷 孔洞 ， 确 认为 缺陷 孔洞 后 
于 逢 闻 周 按 的 特征 绕 园 和 性 测 与 区 本 物 壬 乱 测 区 坊 的 天 生硬 3 之 后 引进 一 个 无 约束 的 三 角 训 分 对 基 曲 面 进行 三 角 
化 ; 最 后 利用 网 格 的 各 向 异性 进行 细 化 及 形态 调整 ， 改 善 网 格 的 拓扑 结构 和 几何 性 质 。 实 验 结果 表明 ， 该 方法 能 够 
有 效 地 识别 三 角 网 格 模型 的 缺陷 孔洞 区 域 并 还 原 其 细节 特征 。 
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Defective hole identification and hole-filling for 3d reconstruction mesh models 


Wen Peizhiab, Lei Yongqinga ?i, Sun Menglonga 
(a. School of Computer Science & Information Security, b. Guangxi Key Laboratory of Intelligent Processing of Images & 
Graphics, Guilin University of Electronic Technology, Guilin Guangxi 541004, China) 


Abstract: Aiming to solve the problem of abnormal holes occur in 3D reconstruction mesh models frequently, this paper 
proposed a novel method that hole-defect identification automatically and detail-preserving surface repair around the hole. 
Firstly, in order to determine whether the hole is a real defective hole, it detected the context irregularity and contour 
irregularity of the missing region. After that, it detected and matched the feature lines around the hole, which were to form 
the fundamental surface of the hole area, and then it introduced unconstrained triangulation for basal surface. Finally, it 
made use of the anisotropy of mesh to refine and adjust morphologically, which were to improve the topological and 
geometric properties of the mesh. The experimental results show that the proposed method can identify the defective hole 
and restore the hole area’s detail features of triangular mesh models effectively. 
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页 河 车 入 马尔 可 夫 随 机 场 (MREF) 进行 规则 化 ; Jug 通 过 建立 一 个 
. 八 叉 树 网 格 结构 ， 提 出 了 等 值 面 曲面 重建 的 方法 ， 但 只 能 适 
在 数字 几何 应 用 领域 ， 三 角 多 边 形 网 格 模型 因 其 良好 的 。 用 处 理 小 的 孔洞 。 这 种 基于 体 的 方法 虽然 在 效率 上 有 良好 的 
可 视 化 、 简 单 灵活 的 表现 能 力 ， 被 越 来 越 多 的 用 在 计算 机 视 。 表现 ,但 是 在 模型 转换 之 间 会 丢失 一 些 原 有 的 细节 特征 信息 ， 
觉 与 图 形 学 中 表示 物体 的 几何 形态 。 由 于 现实 物体 形态 及 材 。 修复 后 的 模型 往往 会 比 原 有 的 模型 出 现 多 余 的 多 边 形 面 片 ， 
质 的 复杂 性 ， 对 目标 物体 进行 三 维 网 格 重建 生成 时 会 造成 孔 ”这 对 后 续 的 网 格 简化 等 操作 会 造成 输出 网 格 质量 的 降低 。 
洞 、 自 相交 、 噪 声 、 边 宛 余 等 方面 的 瑕 羔 ， 这 对 后 续 的 几何 革 于 曲面 定向 的 方法 一 般 是 在 待 修补 的 网 格 模型 上 对 也 
模型 展示 及 应 用 带 来 了 极 大 的 障碍 。 网 格 模型 中 的 孔洞 区 域 。 洞 进行 定位 检测 和 修复 。Liepal9 通 过 最 小 化 相 邻 曲率 的 能 量 
通常 是 模型 缺失 了 几何 信息 且 未 经 前 分 的 缺陷 区 域 ， 待 修补 ，，” 函数 直接 在 待 修复 模型 的 孔洞 区 域 进行 Delaunay 三 角 剖 分 ; 
模型 的 孔洞 区 域 处 没有 任何 已 知 信息 ， 使 得 构造 一 个 连续 、 Wang 等 人 中 先 对 孔洞 顶点 进行 特征 分 类 与 匹配 , 之 后 用 波 前 
合理 的 几何 结构 来 填充 孔洞 区 域 并 得 到 令 人 信服 的 修复 结果 ”法 Advancing Front Method) 进行 细 化 填充 ,， 另 外 还 有 看 
是 一 个 具有 挑战 的 课题 。 根 据 修复 方法 的 不 同 ,Perez 等 人 由 。 于 Bezier 曲面 四 、 径 向 基 函 数 插值 册 、Laplacian 滤波 00、 变 
将 孔洞 修复 的 方法 主要 分 为 两 大 类 : 一 种 是 基于 体 素 的 方法 ， ”分 水 平 集 (1 等 方法 。 这 些 方法 在 一 定 程度 上 能 修复 三 维 模型 
另 一 种 是 基于 曲面 定向 的 方法 。 的 破损 区 域 ， 但 计算 复杂 度 高 且 在 过 渡 区 域 不 能 自然 衔接 ， 
基于 体 素 的 方法 主要 思想 是 将 待 修补 的 网 格 模型 转换 成 。 ”无 法 很 好 的 恢复 孔洞 原始 的 曲面 细节 特征 。 刘 震 等 人 0 将 修 


体 的 形式 表达 ， 在 体 的 基础 上 进行 修补 操作 ， 再 利用 等 值 面 复工 作 分 为 两 步 

由 取 技 术 转 换 表 达成 多 边 形 网 格 。Kumar2 利 用 体 扩 散 方程 面 ， 然 后 在 基 曲 面 的 基础 上 通过 基于 特征 线 优先 级 的 分 层 
治 着 体 的 表面 扩散 到 数据 缺失 区 域 ， 吴 晓 军 等 人 品 提 出 了 本 改 缩 策略 不 断 迭 代 添 加 曲面 细节 ， 直 到 章 分 孔洞 形成 的 三 角 
见 外 壳 与 多 视图 三 维 点 云 融 合 的 多 视图 立体 修复 算法 ; 形 的 边 长 小 于 边界 边 的 边 长 为 止 。 月 ; 
Paulsen 等 人 欠 建 立 符号 距离 场 对 初始 点 集 能 量 最 小 化 ， 后 引 岗 的 每 一 步 算法 复杂 度 高 ， 仪 基于 动态 规划 进行 三 角 剖 分 的 


首先 利用 动态 规划 的 方法 三 角 剖 分 构造 基 
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时 间 复 杂 度 就 为 0(w); 且 现 有 的 大 多 数 方法 中 ,修复 的 过 程 
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验 曲 率 估计 ， 先 查找 当前 点 p(%y,z) 最 邻近 的 ”个 点 ， 构 成 


需要 人 工 参 与 判断 该 孔洞 是 否 是 三 维 重 建 中 形成 的 缺陷 孔洞 
还 是 原本 模型 的 自然 所 属 ， 无 法 满足 自动 化 的 需求 。 
本 文 针 对 上 述 缺 点 ， 提 出 一 种 三 角 网 格 模型 缺陷 孔洞 自 


nx3 的 矩阵 ， 之 后 计算 协 方 差 矩 阵 得 到 特征 值 ， 使 每 一 个 类 
的 方差 都 不 超过 一 个 固定 的 阔 值 65 ， 小 于 阔 值 8 的 部 分 即 为 
选取 的 主 成 分 ， 近 似 于 求解 协 方差 矩阵 进行 主 成 分 分 析 


动 识别 与 曲面 细节 保持 的 修复 方法 ,通过 对 孔洞 的 缺陷 检测 ， (PCA)， 如 下 ; 

解决 了 修复 过 程 中 需要 人 工 判断 的 过 程 ， 而 在 孔洞 填充 部 分 WE 
先 构造 基于 特征 线 的 基 曲 面 ， 之 后 仅 使 用 一 个 时 间 复 杂 度 为 na 

oliogn) 的 三 角 剖 分 函数 ， 不 需 满足 任何 最 优化 的 约 东 条件 ， Is An 

后 续 也 仅 需 较 少 次 数 的 选 代 进行 细 化 ， 避 兔 了 高 代价 的 孔洞 ps 的 名 域 大 小 2 是 每 个 类 的 夺 省 。 风 率 
多 边 形 最 优 三 角 化 的 求解 过 程 。 本 文 方法 中 首先 利用 孔洞 轮 ”估计 如 文献 [13] 计 算 如 下 

廓 周边 的 点 集聚 类 和 轮廓 曲线 切 向 量 的 变化 对 孔洞 进行 检 。 攻 
测 ， 确 认 是 否 为 需要 被 填充 的 单 连通 区 域 缺 陷 孔洞 ， 判 断 为 而 + 

缺陷 孔洞 后 ， 对 缺陷 孔洞 区 域 周边 的 特征 线 进行 检测 与 匹配 。 其 中 ， 为 < 名称， 而 与 入 有 关 的 特征 向 量 表示 数据 的 低 变化 
以 恢复 缺陷 孔洞 区 域 的 特征 线 ， 实 现 缺陷 孔洞 区 域 基 曲 面 的 。 方向 ， 因 此 ， 计 算 表面 p 点 的 法 向 量 ， 便 于 参数 < 能 表示 表 
构造 ， 最 后 对 基 曲 面 填充 时 先进 行 无 约束 的 一 般 三 角 化 ， 再 面 切线 的 变化 ， 这 也 是 表面 区 域 变化 的 一 种 度量 方式 。 只 要 
利用 三 角 网 格 模型 的 各 向 异性 特性 进行 细 化 及 形态 调整 ， 之 。 ”计算 得 到 了 曲线 的 近似 值 ， 那 下 一 步 就 是 计算 表面 的 聚 类 。 


后 能 输出 更 高 质量 、 儿 何 性 质 鲁 棒 的 三 角 网 格 模型 。 实 验 结 
果 表 明 ， 与 文献 [6] 相 比 ， 本 文 方法 很 好 的 保持 了 曲面 原 有 的 
细节 特征 信息 ， 修 复 的 孔洞 区 域 能 与 周边 自然 地 过 渡 衔 接 。 


1 ”缺陷 孔洞 的 检测 


三 维 重建 的 过 程 中 由 于 算法 精度 及 周边 环境 高 光 会 对 重 
建 后 的 三 维 模型 产生 孔洞 的 影响 ， 在 自动 修复 孔洞 过 程 中 需 
要 辨别 该 孔洞 是 需要 被 修复 的 缺陷 孔洞 还 是 模型 自身 所 属 的 
非 缺 陷 孔 洞 。 如 图 1 所 示 , 是 缺陷 孔洞 与 非 缺 陷 孔 洞 的 区 别 。 
自动 修复 孔洞 前 首先 要 搜索 孔洞 的 位 置 ， 而 搜索 孔洞 最 简单 
效 的 方式 就 是 检查 一 系列 边界 边 的 连续 性 ， 所 谓 边 界 边 就 
是 三 角 网 格 模型 中 不 与 别 的 三 角形 共用 同一 条 边 。 文 献 [6] 
中 的 方法 首先 是 在 网 格 上 进行 一 般 搜 索 ， 找 到 一 个 初始 的 边 
界 三 角形 ,接着 对 邻 域 三 角形 进行 搜索 以 寻找 下 一 个 边界 边 ， 
直到 初始 被 检索 ， 一 个 闭环 的 孔洞 路 径 被 检测 出 来 。 
可 是 该 方法 的 局 限 性 在 于 当 被 检测 的 三 维 网 格 模型 是 一 个 非 
完全 封闭 的 表面 时 ， 不 能 判断 此 闭合 路 径 是 否 是 真实 的 缺陷 
孔洞 ， 所 以 自动 修复 孔洞 首先 要 解决 的 问题 在 于 区 分 一 个 表 
面 不 连续 是 否 属于 目标 模型 的 自然 表面 。 基 于 上 述 问题 ， 在 
填充 孔洞 之 前 ， 本 文 首先 判断 网 格 模型 上 孔洞 的 真实 性 ， 提 


以 上 表述 的 算法 流程 如 下 所 示 : 

a) 选 取 一 个 随机 种 子 点 ， 且 生成 一 个 新 的 类 Ci 。 

b) 增 加 一 个 顶点 Ve Nc ， 使 Vo <B 。 这 里 Nc 是 类 ci 的 
径 向 邻 域 ， 且 Vo 是 类 Ci 的 方差 估计 。 

本 文中 用 曲率 估计 的 方法 对 孔洞 轮廓 周边 的 点 集 进行 聚 
类 ， 然 而 高 斯 曲率 和 平均 曲率 这 些 经 典 的 曲面 曲率 表示 方法 
并 不 能 很 好 的 表示 局 部 形态 。 在 这 里 ， 采 用 一 种 独立 于 表面 
大 小 的 测量 方式 ， 以 至 于 能 够 描述 表面 的 局 部 形态 : 形状 索 
引 (shape index)09 就 是 描述 曲面 上 任意 点 在 表面 邻 域 二 级 结 
构 的 一 种 方法 ， 估 计 方 法 如 下 所 示 : 


Ee 


Ss= Fretan fe (3) 
其 中 : kh>k 是 主 曲率 ， 除 了 平坦 区 域 是 一 个 不 确定 的 形状 索 
引 外 ， 其余 表面 区 域 中 se[-b1]。 nm st 
洞 轮廓 ， 就 要 计算 类 间 的 平均 形状 索引 ， 于 形状 索引 
是 一 个 点 方法 ， 所 以 每 一 个 类 i op 

SI= > Dis; (4) 


c 是 K 个 类 中 的 一 个 子 集 ，s 是 类 的 形状 索引 ， 点 集 
的 大 小 N=pi+…+pr ， Pi 表示 每 一 类 的 点 集 数 。 如 图 2 所 示 ， 


出 一 种 验证 缺陷 孔洞 与 非 缺陷 孔洞 的 方法 ， 从 孔洞 最 重要 的 
几何 特性 包括 周边 表面 不 规则 性 和 轮廓 曲线 不 规则 性 这 两 个 
方面 描述 缺陷 孔洞 。 


(a) 缺 陷 孔 洞 
(a)Defect hole 


(b) 非 缺陷 孔洞 
(b)Non-defect hole 
图 1 孔洞 的 分 类 
Fig. 1 Classification of holes 

1.1 孔洞 上 下 文 的 异常 检测 

上 下 文 不 规则 测量 通过 孔洞 周围 区 域 的 差异 量化 进行 聚 
类 环境 的 分 析 。 对 每 一 个 聚 类 来 说 ， 质 心 的 几何 拓扑 特性 用 
来 表示 一 个 形状 特征 ， 因 此 这 一 步 通过 区 域 生长 完成 了 一 般 
性 的 聚 类 生成 。 文 献 [6] 计 算 了 模型 表面 每 一 点 所 在 位 置 的 先 


反映 了 孔洞 周边 点 集聚 类 分 布 的 情况 。 


图 2 ”轮廓 周边 的 点 集 分 布 
Fig.2 Points’ distribution around hole 
1.2 孔洞 轮廓 的 异常 检测 
孔洞 轮廓 的 异常 检测 即 是 检测 计算 轮廓 曲线 的 曲率 与 扭 


力 这 两 个 几何 特性 的 变化 ， 有 具体 参数 的 计算 是 使 
量 的 方法 得 到 。 然 而 基于 


条 件 业 测 
线 扭力 的 方法 ， 不 能 很 好 地 处 理 
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2 ”细节 特征 保持 的 修复 方法 


细节 特 包 


E 的 保持 在 本 文中 指 的 是 恢复 模型 原 有 的 特征 
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当 孔 洞 轮廓 处 于 平坦 区 域 时 的 情况 ， 所 以 在 这 一 步 中 ， 本 文 线 ， 所 谓 三 维 模型 表面 特征 线 是 曲面 上 能 够 表示 显著 几何 特 
采用 文献 [15] 中 对 轮廓 曲线 曲率 测量 的 方法 ， 即 通过 测量 孔 征 的 疹 线 和 谷 线 ， 在 数学 计算 上 的 定义 是 对 模型 表面 沿 着 相 
洞 轮廓 切 向 量 的 变化 ， 判 断 该 孔洞 是 否 属于 缺陷 孔洞 。 对 缺 应 主 方向 曲率 拟 合 极 值 的 描述 。 对 于 具有 丰富 细节 特征 的 
陷 孔 洞 来 说 ， 孔 洞 轮廓 上 切线 的 变化 方向 是 杂乱 无 章 的 ， 反 维 模型 ， 特 征 线 及 其 附近 区 域 经 常 具 有 曲面 的 显著 几何 特征 
映 在 具体 数学 参数 上 表现 为 条 件 和 值 过 高 ;而 非 缺 陷 孔 洞 在 079， 所 以 特征 线 的 检测 与 匹配 连接 对 于 保持 原 孔 洞 区 域 的 细 
切线 角度 上 的 变化 是 平滑 过 渡 的 。 节 信息 具有 重要 的 作用 。 本 文 提出 的 基于 特征 线 的 曲面 细节 
恢复 方法 中 ， 先 进行 特征 线 的 检测 0 与 最 优化 匹配 形成 基 曲 
和 的 构造 ， 这 一 步 完 成 了 孔洞 缺失 部 位 基本 信息 的 恢复 ， 之 
后 对 形成 的 基 曲 面 进行 一 般 的 三 角 剖 分 521， 仅 使 个 算法 
复杂 度 为 O(logn) 的 递归 函数 就 能 实现 , 在 形成 的 基 曲 面 基 础 
上 进行 三 角 面 片 的 痢 分 ， 最 后 利用 孔洞 网 格 的 各 向 异性 进行 
细 化 及 形态 调整 ， 改 善 网 格 的 拓扑 结构 和 几何 性 质 ， 使 之 与 
周边 网 格 完美 过 渡 并 保持 密度 一 致 。 
2.1 特征 线 的 检测 
Sah 特征 线 表征 模型 表面 的 几何 特征 ， 对 带 有 单 连通 区 域 孔 
洞 的 三 维 模型 ， 先 检测 孔洞 周边 区 域 的 特征 线 ， 这 是 本 文 细 
人 @@) 缺 陷 孔 洞 的 法 向 分 布 。。 (b) 非 候 陷 孔洞 的 法 向 分 布 节 特 征 保持 修复 方法 的 第 一 步 。 利 用 离散 微分 算 子 对 分 段 线 
(a)normal distribution of «& (b)normal distribution of non- 性 网 格 进行 操作 ， 避 免 了 高 阶 表面 拟 合 技术 或 隐士 曲面 流程 
defect hole defect hole 这 样 预 处 理 操 作 步 又 。 对 于 绝 大 部 分 几何 模型 来 说 ， 面 噪 
图 3 孔洞 轮廓 的 法 向 分 布 声 造 成 了 大 量 奇异 三 角形 的 出 现 ,在 进行 特征 线 的 检测 之 前 ， 
Fig.3 ”Normal distribution of hole boundary 可 先 对 三 维 表面 进行 平滑 处 理 操作 ， 之 后 再 计算 离散 极 值 。 
切 向 量 定 义 为 TO=NGSxBSsS=rOlrl ， 其 中 假设 M 是 三 维 空 间 中 的 光滑 曲面 ，K; 是 沿 着 M 中 曲率 
N(s)=r"(D/r"O 为 法 向 量 ， sd = ACO 为 参数 曲线 "上 的 线 上 的 单位 切 向 量 , @ =(VKi,K) 被 称 为 是 极 值 系 数 。 其 中 e 是 
2 ey Re i Ub PN Cy 曲面 上 的 隐 函 数 ，K 的 符号 选择 决定 了 的 符号 。 对 于 M 上 
i se ge pl he 
采用 最 小 二 乘法 和 局 部 弧 长 估计 的 方法 来 计算 其 切 向 量 ， 根 数 e 的 值 为 0 即 为 曲面 M 上 的 一 条 特征 线 ,分 为 两 种 情况 em 
据 条 件 炉 9 的 概念 ， 切 线 变 化 的 条 件 业 为 和 em : 
Si =-》 po)log(p(o (5) evn = 0: (Vewws Rann) < OK | > IK 0) 
综 上 所 述 ， 孔 洞 轮廓 异常 测量 (全 局 粹 ) 的 计算 方式 为 I 加 
irregularity =||SI|: Ss (6) 网 二 、 、 
对 于 规则 三 角形 T=(p,ps,p;) ,依次 对 三 角形 的 顶点 计算 
Kwa(pi) 和 emw(p) ，i=(1,2,3) 。 由 于 ew 对 于 三 角形 T 是 一 个 线 
性 函数 ， 计 算 e。 的 梯度 。 
B (Vess Daer Ka (pi)) <0 (9) 
[Er Ros (p)]> | Rs (10) 
x 如 果 eu =0 为 非 空 且 满足 式 (9)(10) 的 条 件 , 与 ee 对 应 的 
N 特征 线段 就 被 检测 出 来 。 对 式 (8) 采 取 上 述 类 似 的 处 理 方式 ， 
则 对 应 em 的 特征 线段 也 被 检测 出 来 。 在 所 有 的 特征 线 被 提 
取 之 后 ， 用 闵 值 方案 处 理 所 提取 的 特征 线 ， 按 弃 一 些小 线段 
和 弯曲 的 波浪 线 ， 余 下 的 即 为 本 文 方法 所 采用 的 特征 线 。 
2.2 ”特征 线 的 匹配 
图 4 曲线 上 切线 与 法 线 的 分 布 对 三 维 网 格 模 型 孔洞 周边 的 特征 线 进行 了 检测 ， 接 下 来 
Fig.4 Tangent, normal distribution of curve 就 是 进行 特征 线 的 匹配 2 以 恢复 孔洞 中 的 特征 线 ， 构 造 三 角 
图 3 所 示 是 缺陷 孔洞 与 非 缺 陷 孔 洞 关 于 切 向 量 方向 的 变 剖 分 的 基 ， 以 提高 孔洞 区 域 上 的 细节 恢复 效果 。 
化 情况 ， 图 4 所 示 是 曲线 上 切 向 量 与 法 向 量 的 分 布 。 实 验 中 假设 5={s},i=1,2,3,…N 是 孔洞 区 域 周边 提取 的 特征 线 集 
设 对 参照 实验 各 10 个 孔洞 轮廓 , 一 组 是 缺陷 孔洞 ， 人 合 ， 对 于 特征 线 sss ， 待 匹配 特征 线 集 合 为 
组 是 非 缺陷 孔洞 。 经 计算 得 到 使 缺陷 孔洞 与 非 缺陷 孔洞 分 本 
的 确信 确定 为 3 ， 调 王 3. 的 病人 确定 为 全 了 孔 油 ， 反 之 为 全 55 局内 菱帅 5 是 区 和 条 过 乔 训 训 吕 尖 起 吉 有 有 
非 缺 陷 孔 洞 。 


vw 的 特征 线 ， 引 = 这 (wa -vin) 为 特征 线 


和 位 方向 ， 阀 值 参数 3e[-V2/2,0] ， 匹 配 集合 中 元 素 的 数 


录用 定稿 
目 决定 于 阀 值 5 的 大 小 ，5 越 小 则 集合 中 元 素 的 数目 就 越 多， 
实际 中 本 文选 取 5 的 值 为 -V2/2 。 给 出 了 特征 线 5 与 集合 c 中 


的 待 匹 配 特征 线 ， 通 过 匹配 概率 (matching possibility,，MP) 
衡量 标准 选取 最 优 匹配 : 


] 才 Wo 
pees) -| | (11) 


CO2 


其 中 : 5s;eG，an 为 5; 的 第 k 个 采样 点 vi 在 特征 线 s 拟 合 曲 
线 上 的 投影 , 高 斯 函数 中 参数 o 取 0.5, 如 图 5 所 示 为 特征 线 
检测 与 匹配 的 示意 图 。 

对 于 特征 线 s， 使 MP(s,s)) 最 大 的 5 sc 即 为 s 的 最 优 匹 
配 。 如 果 s 的 最 优 匹 配 也 是 s ， 则 称 (5,,5)) 为 最 优 匹配 对 ， 选 
取 最 优 匹 配对 (5.,5)) ,使 用 球面 线性 插值 方法 将 特征 线 匹配 对 


OP (12) 


其 中 : ,=(Vio —0), 于 =(vjo -0)， 其 中 o 是 


确定 圆 的 圆心 ，9 为 向 量 布 和 避 之 间 的 夹 角 ， 连 接 起 来 的 这 
条 线 称 为 匹配 对 (si,5)) 的 桥 线 。 


ViosVisvVio 三 点 所 


图 5 特征 线 的 检测 与 匹配 
Fig.5 Detection and matching of feature lines 

2.3 三 角 训 分 及 细 化 调整 

上 述 部 分 进行 特征 线 的 匹配 与 桥接 之 后 ， 表 征 细节 特 
信息 的 桨 线 与 谷 线 被 恢复 ， 原 孔洞 部 分 被 划分 为 了 几 个 简 
的 多 边 形 子 洞 ， 基 曲面 构造 完成 。 由 于 形成 的 基 
征 线 的 连接 ， 并 没有 形成 表示 网 格 模 型 的 三 角 面 片 ， 之 后 需 
对 每 一 个 多 边 形 子 洞 进 行 三 角 化 。 现 有 的 许多 方法 中 ， 和 孔洞 
区 域 的 三 角 齐 分 是 整个 算法 的 关键 步 又， 其 本 质 是 一 个 非 线 
性 的 优化 问题 ， 为 获得 最 优 的 三 角 化 结果 ， 其 时 间 复 杂 度 通 
常 很 高 。 为 了 解决 最 优化 三 角 训 分 时 间 复 杂 度 高 的 问题 ， 在 
本 方法 之 前 步骤 中 形成 的 基 曲 面 的 基础 上 使 用 平均 时 间 复 杂 
度 为 ol(logn) 的 递归 函数 Find(i,j) 进行 无 约束 的 三 角 化 (i ，j 
为 多 边 形 子 洞 的 项 点 标号 )， 步 又 所 示 如 下 : 

a) 如 果 j=i+2 ， 则 将 Avivinvj 加 入 到 三 角形 集合 s 中 。 

b) 否 则 令 mi 为 i，j 之 间 的 任 一 数 。 

c) 在 步骤 b) 的 条 件 下 ， 如 果 mi 不 等 于 i+1， 递 归 调 用 


Find(i,my) ， 将 Awywv 加 入 到 三 角形 集合 s 中 。 
中 如 果 m 不 等 于 j-1， 弟 归 调 用 FindCm，, 旋 
在 递归 函数 Find(i, 站 中 ，Avivwvj 是 本 次 递归 分 割 出 的 三 
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晶 此 时 填充 的 网 格 质量 较 差 ， 与 周围 网 格 的 采样 率 及 顶点 密 
度 相 差 较 大 ， 所 以 需要 对 填充 后 的 孔洞 网 格 进行 细 化 。 本 文 
采用 边 记分 函数 为 每 条 边 赋予 一 个 影响 因子 09， 实 现 网 格 的 
各 向 异性 的 细 化 。 

影响 因子 由 边 的 长 度 和 密度 两 个 元 素 构成 ， 反 映 的 是 最 
终 的 孔洞 网 格 的 疏 密 程度 ， 边 ai=wrw， 则 边 的 长 度 为 


l(ej)=|les|=|v, -vl 。 对 于 男 一 个 元 素 边 的 密度 ， 先 定义 顶点 v 


的 密度 pw =wKCI+d-AH2) ，ye(0,D) ， 而 K(v) 和 瑟 (V) 分 
别 为 该 点 的 高 斯 曲率 和 平均 曲率 。 如 下 图 6 所 示 ， 对 于 顶点 
Ww ， enen en 为 其 顶点 w 对 相 邻 顶点 vw…,v 形成 的 边 ; 
相 邻 两 条 边 e 和 eo 之 间 的 夹 角 为 0;=ang(ew,e0w) ; 


Eo0i X eo(inl 村 
Avovoivouin 所 在 平面 的 法 向 量 为 "= 相 邻 三 角形 所 


lev: xeooan| 
在 平面 的 夹 角 为 Pi =ang(ni,nin) 。 
顶点 w 的 高 斯 曲率 K 和 平均 曲率 号 BO 的 定义 为 


~ 


a 


2r -Ya 
K(W)=— Q3) 
3area(S) 
>lellal 
Ho)= 气 -一 (14) 
4area(S) 
其 中 :Ss 是 与 顶点 内 相 邻 的 所 有 三 角形 的 面积 之 和 。 


OW 


图 6 与 顶点 w 相关 的 几何 量 定义 
Fig.6 Geometric definition about vo 
本 文 假设 边 的 密度 从 一 个 顶点 v 到 另 一 个 顶点 vj 之 间 
是 线性 变化 , 则 边 密 度 clej) 可 由 两 个 顶点 密度 的 均值 表示 为 
Ea) (15) 


Ple;)= 


定义 了 边 的 长 度 和 密度 之 后 ， 边 @ 的 记分 函数 由 准 高 其 

核 函 数 给 出 ， 即 
| 

w= p73 )] (16) 

通过 该 记分 函数 ， 先 计算 孔洞 周围 区 域 网 格 的 每 条 边 的 

记分 ， 并 求 其 平均 值 页 。 和 孔洞 网 格 填充 后 的 每 条 边 的 一 个 记 

， 若 My> 斩 ， 则 对 该 边 进行 边 分 裂 操 作 ， 即 在 该 边 

的 中 点 处 插入 一 个 新 的 项 点， 并 将 这 个 新 顶点 与 该 边 所 在 三 


角形 的 两 个 对 顶点 分 别 相 连 ， 如 图 7 所 示 。 
对 孔洞 网 格 填 充 细 化 之 后 ， 需 要 对 修复 后 的 网 格 进行 几 


角形 ， 依 次 递归 ， 多 个 多 边 形 子 洞 被 三 角 剖 分 而 得 以 填充 。 


何 形态 调整 ， 这 种 调整 在 局 部 是 按 各 向 异性 进行 的 。 本 文通 
过 边 置 换 来 保证 孔洞 网 格 的 Delaunay 性 质 。 对 于 相 邻 的 两 个 
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三 角形 AABC 和 AABD ， 当 


有 局 部 Delaunay 性 质 ， 无 须 进行 操作 ;， 和 否则 需要 将 边 4B 置 
换 为 边 CD ， 如 


熙 己 
Ti 


3 


算法 的 有 效 性 与 正确 性 ， 实 验 对 象 选择 图 形 库 Geometry++ 
尺 表 的 三 维 网 格 结构 ， 修 复 后 的 结果 也 均 以 网 格 模型 
展示 ， 实 验 结果 与 Liepa 方法 进行 对 比 。 


中 取 
的 结构 


图 8 所 示 。 


AABC 确定 一 个 外 接 圆 时 , 对 顶点 D 
进行 判断 : D 到 该 圆 上 的 投影 点 D' 在 外 接 圆 外 时 ， 边 AB 具 


-< 


图 7 边 的 拓扑 分 裂 
Fig.7 Edge’s topology split 


Fig. 8 


图 8 边 的 拓扑 置换 

Edge’s topology substitution 
上 述 过 程 在 孔洞 区 域 一 般 三 角 化 后 进行 迭代 操作 ， 直 到 

无 边 分 裂 即 完成 了 孔洞 网 格 的 细 化 。 这 种 后 续 的 规整 步骤 只 


分 别 为 图 (a)(b)(c) 孔 洞 
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图 9 所 示 为 Eight 模型 的 修复 效果 对 比 ， 图 中 ，(d)(e)(f) 
右 旋转 60 度 后 的 观 


察 角 度 ，(a)(d) 为 原始 破损 模型 ，(b)(e) 为 Liepa 方法 的 结果 ， 


(上 为 本 文 方法 修复 后 的 效果 图 。 
方法 只 是 在 孔洞 区 域 原始 顶点 的 基础 上 进 


从 图 


中 可 以 看 出 ，Liepa 


行 三 角 训 分， 破坏 


了 原始 模型 的 几何 结构 ;而 本 文 方法 恢复 了 孔洞 区 域 的 几何 


特征 结构 并 且 延 续 了 孔洞 周边 区 
线 特征 得 到 了 很 好 的 保持 。 图 
图 中 ，(d)(e)( 分 另 


结果 对 比 ， 


成 的 特征 线 ， 对 其 模型 的 兰 
10 所 示 为 Squirrel 模型 的 修复 
上 为 图 (a)(b)(c) 的 局 部 放大 图 ， 


(b)(e) 为 Liepa 方法 的 结果 图 。 


格 密度 与 周边 三 角 网 格 密度 不 一 致 且 没有 


可 以 明显 看 出 修复 后 的 三 角 网 


民 好 的 延续 原 有 了 筷 


洞 周 边 网 格 的 特征 线 ，(c)(f) 为 本 文 算法 修复 后 的 效果 ， 不 仅 


修复 后 的 三 角 网 格 密度 与 
致 ， 而 且 很 好 的 延续 了 孔洞 周边 的 特征 线 ， 


原始 孔洞 周边 三 角 网 格 密度 保持 一 


获 了 孔洞 的 全 


局 特征 ， 使 孔洞 区 域 的 显著 性 几何 特征 得 到 了 很 好 的 恢复 。 


11 所 示 为 Diana 模型 的 修复 效果 比较 , (a) 为 带 孔 洞 的 原始 
模型 ， 该 模型 的 后 脑 顶 部 有 大 部 分 缺失 区 域 ，Liepa 方法 只 


能 沿 着 孔洞 轮廓 项 点 进 


后 脑 顶部 凸 起 特 和 


行 三 角 放 分， 根本 不 能 还 原 原 始 模型 
E， 修 复 效果 极 差 ， 难 以 反映 其 原始 模型 的 


原始 特征 信息 ; 本 文 方法 修复 后 的 效果 可 以 明显 看 出 对 后 脑 


顶部 的 孔洞 区 域 原始 特征 
密度 也 与 原始 网 格 模型 保持 一 致 。! 
出 ， 本 文 方法 的 修复 结果 不 仅 能 捕获 孔洞 区 域 的 全 局 特征 ， 


言 息 进行 了 很 好 的 还 原 ， 三 角 网 格 


要 较 少 次 数 的 简单 迭代 ， 就 能 获得 几何 怕 
格 ， 且 能 够 与 孔洞 周边 实现 自然 过 渡 。 


实验 结果 与 分 析 


FE 质 良好 的 三 角 网 


本 文 实验 电脑 配置 为 : CPU 3.70 GHz Intel Core™ 
i3-6100、RAM 12.0 GB、Windows10 64 位 ， 编 程 采用 
MATLAB 、Python，OpenCV、MeshLab 。 为 验证 本 文 提 昌 


的 


且 局 部 的 显著 几何 特 和 
孔洞 区 域 都 给 出 了 较 好 的 填充 效果 。 图 


了 | 


下 也 得 到 了 很 好 


实验 的 修复 效果 可 以 看 


的 恢复 ， 对 不 同 尺 度 的 
12 中 为 Bunny 模型 


并 24 


的 修复 效果 比较 ，(a) 为 带 两 个 单 连通 区 域 的 孔洞 模型 ， 其 中 


大 腿 和 背部 都 有 明显 的 区 域 


复 效 果 跟 图 


11 所 示 相 


决 失 ; 经 过 Liepa 方法 后 ， 其 修 


司 , 只 能 沿 着 孔洞 轮廓 顶点 进行 三 角 齐 


分 ， 且 没有 考虑 到 采样 密度 与 周边 三 角 网 格 一 致 ， 修 复 效 果 


极 差 , 难以 反映 模型 


且 在 Bunny 背部 及 
细节 特征 保持 良好 。 


原 表 面 的 细节 特征 信息 , 如 图 
图 12(c) 为 本 文 的 修复 方法 , 不仅 还 原 了 颈 部 凹陷 的 特征 ， 并 
其 大 腿 上 很 好 的 修复 了 原 有 的 凸 起 特征 ， 


12(b) 所 示 ; 


(a)Eight 模 型 


(a) The original Fight model 


(da 


图 的 局 部 放大 


(b)Liepa 方 法 
(b)Repaired by Liepa 


闻 
(©)b 
(d) The local enlarged drawing of a (e) The local enlarged drawing of b 
图 9 Eight 模型 修复 效果 比较 
Fig.9 Comparison of Eight model 


图 的 局 部 放大 


(0) 本 文 方法 


(c)Repaired by this paper 


(fc 图 的 局 部 放大 
(f) The local enlarged drawing of c 
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(a)Squirrel 模 型 (b)Liepa 方 法 (c) 本 文 方法 
(a) The original Squirrel model (b)Repaired by Liepa (c)Repaired by this paper 
有 


(da 图 的 局 部 放大 (e)b (Dc 图 的 局 部 放大 
(d) The local enlarged drawing of a (e) The local enlarged drawing of b (f) The local enlarged drawing of c 


图 10 ”Squirrel 模型 修复 效果 比较 


Fig. 10 Comparison of Squirrel model 


(a)Diana 模 型 (b)Liepa 方 法 (c) 本 文 方法 
(a) The original Diana model (b)Repaired by Liepa (c)Repaired by this paper 


图 11 Diana 模型 修复 效果 比较 


Fig. 11 Comparison of Diana model 


计 


(a)Bunny 模 型 (b)Liepa 方 法 (c) 本 文 方法 
(a) The original Bunny model (b)Repaired by Liepa (c)Repaired by this paper 


图 12 Bunny 模型 修复 效果 比较 


Fig. 12 Comparison of Bunny model 
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表 1 和 2 分 别 列 出 了 原始 模型 、Liepa 方法 和 本 文 方法 
所 修复 后 的 模型 顶点 个 数 与 三 角 面 片 的 数量 比较 ， 
Liepa 方法 后 的 顶点 数量 与 原始 破损 模型 的 顶点 数量 一 样 ， 
即 Liepa 方法 只 在 原始 破损 模型 的 孔洞 轮廓 上 进行 三 角 齐 6 


2010, 16 (4): 636-646. 


Ju Tao. Robust repair of polygonal models [J].ACM Trans on Graphics. 
2004, 23 (3): 888-895. 
Liepa P. Filling holes in meshes [C]/Proc of Eurographics/ACM 


分 ， 不 会 考虑 到 原始 特征 信息 的 保持 ， 而 本 文 方法 修复 后 以 SIGGRAPH Symposium on Geometry Processing. 2003: 200-205. 
模型 不 管 是 在 顶点 数量 还 是 在 三 相片 的 数量 上 都 高 于 7] Wang Xiaochao, Liu Xiuping, Lu Linfa, et al. Automatic hole-filling of 
以 人 3 1 可 8 ping 8 
Liepa 方法 ， 且 对 于 大 孔洞 区 域 的 面 片 数 大 大 增加 ， 极 大 的 CAD models with feature-preserving [J]. Computers & Graphics, 2012, 
丰富 了 孔洞 区 域 的 细节 特征 信息 。 36 (2): 101-110. 
表 1 模型 的 顶点 数 比较 8] Li Zhong, Meek D S$, Walton D J. Polynomial blending in a mesh 
Table ] Vertexes comparison of models hole-filling application [J]. Computer-Aided Design, 2010, 42(4): 
原始 模型 Liepa 方法 本 文 方法 340-349. 
Eight 766 766 768 9] Liu Shengjun， Wang C C L. Quasi-interpolation for surface 
Squirrel 10745 10745 10850 reconstruction from scattered data with radial basis function []]. 
Diana 20 910% L004 Computer Aided Geometric Design, 2012, 29(7): 435-447. 
Bunny 33144 33144 35292 ye a 证 A 6 
二 开 | 中 六 2 el ingqlang, u Jianhuang， 上 an ingyong. n integrate' 
表 2 模型 的 面 片 数 比较 : pa ee . 
. approach to filling holes in meshes [Cl]// Proc of International 
Table 2 Patches comparison of models 人 i i 1 
二 a i Conference on Artificial Intelligence and Computational Intelligence. 
原始 模型 ” ”Liepa 方 法 本 文 方法 8 
Eight 1516 1532 1543 ee 
Squirrel 21454 21465 21477 [11] 张 婧 ， 周明 全 , 耿 国 华 .基于 变 分 水 平 集 的 三 维 模型 复杂 孔洞 修复 
Diana 19379 19520 20074 []. 计算 机 应 用 研究 ，2018，35(4):1254-1257.(Zhang Jing，Zhou 
Bunny 65922 66280 70580 Mingquan, Geng Guohua. Repair methods for complex hole of 


three-dimensional model based on variational level set method [J]. 
Application Research of Computers, 2018, 35(4): 1254-1257. ) 

[12] 刘 震 , 王 艳 宾 ， 白 丽 丽 ,， 等 . 曲面 细节 特征 保持 的 三 维 模型 孔洞 修复 
方法 [J]. 计算 机 辅助 设计 与 图 形 学 学 报 , 2016, 28(12): 2052-2059. 


4 ”结束 语 


针对 任意 三 角 网 格 模 型 ， 本 文 提出 了 一 种 自动 识别 真实 
缺陷 孔洞 与 保持 细节 特征 的 修复 方法 ， 解 决 了 现 阶 段 大 多 数 
孔洞 填充 算法 都 需要 人 工 参与 判断 的 问题 ， 利 用 无 约束 的 三 
角 化 代替 求解 非 线性 的 最 优化 三 角 放 分 极 大 地 降低 了 算法 的 
时 间 复 杂 度 。 首 先 利 用 孔洞 上 下 文 异 常 检测 及 轮廓 异常 检测 
对 三 角 网 格 模型 的 孔洞 进行 检测 ， 判 断 该 缺失 区 域 是 否 属于 13 


真实 缺陷 孔洞 ， 然 后 利用 孔洞 周边 特征 线 的 信息 进行 检测 与 


T 


(Liu Zhen, Wang Yanbin, Bai Lili, et al. Detail-preserving hole-filling 


for complex 3D models [J]]. Journal of Computer-Aided Design & 
Computer Graphics, 2016, 28(12): 2052-2059. ) 
Pauly M, Gross M, Kobbelt L P. Efficient simplification of 


point-sampled surfaces [Cl]//Proc of Conference on Visualization. 


匹配 以 实现 对 破损 区 域 基 曲 面 的 构造 ， 再 在 基 曲 面 的 基础 上 Washington DC:IEEE Computer Society, 2002: 163-170. 

进行 三 角 痢 分 ,采用 边 记分 函数 为 每 条 边 赋予 一 个 影响 因子 ， 14] Koenderink J J, Van Doorn A J. Surface shape and curvature scales [J]. 
进而 实现 网 格 的 各 向 异性 的 细 化 及 形态 调整 。 通 过 与 其 他 方 Image and vision computing, 1992, 10(8): 557-564. 

法 相 比 ， 本 文 能 很 好 地 还 原 原 始 模型 的 旨 节 特 和 上 且 在 模 15] Sanchez G, Branch J W, Atencio P. A metric for automatic hole 


型 修复 质量 和 计算 效率 方面 取得 了 很 好 的 平衡 ， 为 三 维 模型 
的 后 续 处 理 提供 更 精确 的 模型 结构 。 但 本 文 解决 的 问题 仅 适 
; 失 
步 


characterization [Cl]//Proc of the 19th International Meshing Roundtable. 
Berlin:Springer,2010: 195-208. 


用 于 单 连通 区 域 的 缺陷 孔洞 修复 ， 忽 略 了 复杂 缺失 区 域 的 拓 16 


Caticha A. Lectures on probability, entropy, and statistical physics 


扑 结构 ， 算 法 的 时 间 复 杂 度 也 较 高 ， 后 期 会 进 研究 复杂 [Cl/Proc of International Workshop on Bayesian Inference & 

缺失 区 域 的 情况 ， 考虑 使 用 GPU 并 行 加 速 处 理 ， 以 实现 复 Maximum Entropy Methods in Science & Engineering.2008. 

杂 破 损 模型 的 快速 修复 。 17] Ohtake Y, Belyaev A, Seidel H P. Ridge-valley lines on meshes via 
ks implicit surface fitting [J].ACM Trans on Graphics,2004，23(3): 
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